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Abstract. We consider a minimum-weight perfect matching problem on the 
line and establish a "bottom-up" recursion relation for weights of partial 
minimum-weight matchings. 



1. Introduction 

We start with recalling a few notions from combinatorial optimization on graphs. 
A matching in an undirected graph is any set of its mutually disjoint edges: no two 
edges from such set can share a vertex. A matching is called perfect if it involves 
all vertices of the graph (the number of vertices is then necessarily even). 

Depending on the structure of the graph, perfect matchings may be many. Sup- 
pose that edges of a graph are endowed with real weights] then it makes sense to 
look for a perfect matching composed from a set of edges with a minimum sum 
of weights. In this note we treat a particular case of this minimum-weight perfect 
matching problem where the graph is complete, all its vertices are located on a line, 
and edge weights are related to distances between the vertices along the line. 

A bipartite version of this problem, in which vertices are divided into two equal 
classes and edges of the matching must connect vertices of different class, reduces to 
transport optimization. For a particular class of cost functions (those of "concave 
type") this problem has been thoroughly treated in the measure-theoretic setting 
in |10j . Similar problems have also long been considered in the algorithmics liter- 
ature for the specific case of the distance \x — y\, assuming that the measures are 
discrete [Uinilll]- An algorithm for a general cost function of a concave type has 
been proposed recently in [31[S] . 

In the present note we consider the minimum-weight perfect matching problem 
on a complete graph without assuming it to be bipartite. However, for the class 
of weight functions generated by distances, the non-bipartite problem turns out to 
be essentially equivalent to a bipartite problem with alternating points due to a 
"no-crossing" property of the optimal matching. This observation, which we owe 
to S. Nechaev, allows to employ the theory developed for the bipartite case in [4l[5] 
to the non-bipartite matching. 

The main contribution of the present note is a specific "bottom-up" recursion 
relation ([6|) for partial minimum- weight matchings. This recursion follows from a 
"localization" property of minimum-cost perfect matchings for concave cost func- 
tions (Theorem [5]) and provides a new perspective on the construction of [HH]. 
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The note is organized as follows. Section [2] is a review of the basic construc- 
tion of metrics on the real line. In Section [3] we recall the problem of minimum- 
weight perfect matching and cite a few useful results about the structure of its 
solution. Section |4] contains the key technical result, a kind of localization prin- 
ciple for minimum- weight matching. Finally in Section [5] we derive the recursive 
relation for partial minimum-weight matchings. We also discuss the ensuing algo- 
rithm and compare it to modern variants of the Edmonds blossom algorithm for 
the minimum-weight perfect matching problem. 

It is our pleasure to thank the organizers of the conference Optimization and 
Stochastic Methods for Spatially Distributed Information (St Petersburg, EIMI, 
May 2010), where an earlier version of this work was presented. 



Recall the usual axioms for a distance ■) on the real line R: for all x,y, z G R, 

(Dl) y)>0 with d{x, y) ^ iS x ^ y; 

(D2) d{x,y) = diy,x); 

(D3) d{x,y) + d{y, z) > d{x, z). 

These axioms are satisfied by the distance d{x, y) = |a; — ?/| as well as by any distance 
of the form 

(1) dg{x,y) = g{\x-y\), 

where g{-) is a nonnegative concave function defined for all a; > such that g{x) — 
iff a; = 0. Here concavity means that g{Xx -I- (1 — X)y) > Xg{x) -f (1 — X)g{y) for all 
x,y >0 and all A G [0, 1]. Note that the distance dg is homogeneous with respect 
to shifts: 

(D4) dg{x + t,y + t) = dg{x, y) for all x,y,t G R. 

Conversely, any distance satisfying the axioms (E[T])-(E|4]) has the form ([T|) for a 
suitable nonnegative concave function g: indeed, take g{x) — d{0,x) for a; > and 
check that concavity follows from (E(3]). 

An example of this construction is the distance \x — yl" with < a < 1, where 



is the "discrete distance"; here we are mostly interested in the case < a < 1. 

An important property of a distance d{-, ■) is whether it is intrinsic. To recall 
the corresponding definition, take two distinct points a;, y S R and connect them 
with parameterized curves taking values x{t) in R, i.e., suppose that < t < 1, 
x(0) = X, and x{l) = y. By the triangle inequality, 



where the infimum is taken over all curves connecting x to y and all meshes = 
to < ti < t2 < ■ ■ ■ < t]^ = 1 with TV > 1. The distance d is called intrinsic if ([2]) is 
an equality for all x,y (see, e.g., [2]). 

The distance |a; — j/| and its scalar multiples are the only homogeneous intrinsic 
distances in R. For distances \x — y\" with < a < 1, connecting curves have 
infinite length unless x — y, and therefore these distances are not intrinsic. 



2. Intrinsic and non-intrinsic distances on the line 




(2) 




0<i<JV 
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The intuition behind this notion is that whereas the geometry of R equipped 
with an intrinsic distance is fuhy determined by the distance, the same hne R with 
a non-intrinsic distance should be viewed as embedded in an auxihary space of a 
larger dimension, and geometry on R is induced by that in the embedding space. 
According to the Assouad embedding theorem (see, e.g., [8]), for the distance \x—y\°' 
on R the embedding space has dimension of the order 1/a for small a. 



Consider an even number of points xi < X2 < ■ ■ ■ < X2n on the real line R 
equipped with a distance d and look for a minimum-weight perfect matching in a 
complete graph i^2n on these points, where the weight of an edge connecting Xi 
to Xj is d{xi ,Xj). It is convenient to represent R with a horizontal interval and use 
arcs in the upper halfplane to show the edges of the graph i^2n- 

The following lemma, adapted from McCann [TUl Lemma 2.1], allows to describe 
the structure of minimum- weight matchings. 

Lemma 1. Suppose the distance has the form dg ([1]) with a strictly concave func- 
tion g (i.e., g{Xx + (1 — X)y) > Xg{x) + (1 — X)g{y) for all x,y > 0, x y, and 
< X < 1). Then the inequality 

(3) dg{xi,yi) -\- dg{x2,y2) < dg{xi,y2) -\- dg{x2,yi) 

implies that the intervals connecting Xi to yi and X2 to y2 are either disjoint or 
one of them is contained in the other. 

Proof. Let arcs representing the matching be directed from x's to j/'s and assume, 
without loss of generality, that xi is leftmost. Then the configuration 
where xi < y2 < yi < X2, cannot satisfy (jS]) because a strictly concave function g 
must be strictly growing and therefore dg{xi,y2) < dg{xi,yi) and dg{x2,yi) < 
dg{x2, 2/2)- To rule out configuration where xi < X2 < yi < y2, choose 

< A < 1 such that 

yi - xi = (1 - X){y2 - xi) + X{yi - ^2). 
Then, since (yi - xi) + (2/2 - X2) = {y2 - xi) + [yi - X2), we obtain 

y2- X2 = X{y2 - xi) + (1 - A)(?;i - X2) 
and can further use the strict concavity of g to get 



dg{xi,yi) > (1 - X)dg{xi,y2) + Xdg{x2,yi), 

dg{x2,y2) > Xdg{xi,y2) + (1 - X)dg{x2,yi)- 

The sum of these inequalities contradicts ([3]). All the other configurations where 
the points X2, yi, y2 are located to the right of xi, namely 



It is possible to show that conversely, for an arbitrary bivariate function d(x,y) 
the property of Lemma [T] together with homogeneity (E|3I) imply that d has the 
form ^ for a suitable concave function g [TO] Lemma B4] . 

We call a matching nested if, for any two arcs {xi,Xj) and {xii,Xj>) that are 
present in the matching, the corresponding intervals in R are either disjoint or one 
of them is contained in the other. 

Theorem 2 (plllOj). A minimum- weight matching is nested. 



3. Minimum-weight perfect matchings 




□ 
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Proof. This statement essentially follows from Lemma [TJ because whenever two 
arcs {xi,Xj) and {xi',Xj') are crossed, the sum of distances corresponding to the 
"uncrossed" arcs {xi,Xji), {xi',Xj) must be smaller according to ([3]). 

However uncrossing may introduce new crossings with other arcs, and it remains 
to be checked that it does decrease the total number of crossings. This is done 
using an argument adapted from l,, Lemma 1]. 

Let an arc {xk,xi) cross either of the arcs {xi,Xj>) or {xi',Xj) obtained after un- 
crossing the arcs {xi, Xj) and {xi' , Xj'); we would like to prove that before uncrossing 
the arc {xk,xe) had the same number of crossings with (xi,Xj) and {xi',Xj'). It 
suffices to consider the following two cases. 

Case Xi < Xj' < Xj < Xi'. Assume that {xk,xe) crosses {xi,Xj'). Then {xk,X(), 
shown as a dotted arc, may be situated in one of the following ways: 



after uncrossing before uncrossing 




tXj ^ J 3 ^ ^ 3 3 ^ 



We see that the number of crossings between {xk,Xi) and the other two arcs is 
the same before and after uncrossing. The case where {xk,Xi) crosses {xi',Xj) is 
symmetrical. 

Case Xi < Xii < Xj < Xji. The following situations of the dotted arc {xk,Xi) 
are possible: 



after uncrossing before uncrossing 




Again the number of crossings between (xk, xe) and the other two arcs is the same 
before and after uncrossing. 

It follows that each uncrossing removes exactly one crossing from the matching, 
and therefore any possible sequence of uncrossings leads in a finite number of steps 
to a nested matching with a strictly smaller weight. In other words, it suffices to 
look for the minimum-weight matching only among the nested ones. □ 

This result implies that the minimum-weight perfect matching problem is es- 
sentially bipartite, and it is indeed the bipartite setting that is considered in 
[Illl[5l[9Hll] 

Corollary 3. In a minimum-weight perfect matching, points with even numbers 
are matched to points with odd numbers. 

Proof. Suppose on the contrary that Xi is matched to Xj with i and j both even or 
both odd. The interval connecting Xi to Xj contains an odd number |i — j| — 1 of 
points, one of which has to be matched to a point outside and thus cause a crossing 
with (xi,Xj). Since a minimum- weight matching is nested, the result follows. □ 

Observe finally the following simple form of the Bellman optimality principle. 
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Vk ■■■ yi XI zx ■■■ zk z'f., ■■■ z[ Xr y'l ■■■ y'k' 



Figure 1. Notation used in the proof of Theorem [5] Note that in 
general Xi < zi and z[ < Xr, but in the figure these pairs of points 
are shown to be distinct. 

Lemma 4. Any subset of arcs in a minimum-weight matching is itself a minimum- 
weight perfect matching on the set of endpoints of the arcs that belong to this subset. 

Proof. Indeed, if one could rematch these points achieving a smaller total weight, 
then the full original matching itself would not be minimum- weight: rematching the 
corresponding subset of arcs (and possibly uncrossing any crossed arcs that might 
result from the rematch) would give a matching with a strictly smaller weight. □ 

4. Preservation of hidden arcs 

Call an arc {xi,Xj) in a nested matching exposed if it is not contained in any 
other arc, i.e., if there is no arc {xii,Xj>) with Xi,Xj contained between Xi' and Xji. 
We call all other arcs in a nested matching non-exposed or hidden. Intuitively, 
exposed arcs are those visible "from above" and hidden arcs are those covered with 
exposed ones. 

Suppose X = {Xi}i<i<2n with Xi < X2 < ■ ■ ■ < X2n and X' = {x^,}i<t'<2n' 

with x[ < x'2 < ■ ■ ■ < x'^^i be two sets such that X2n < x[, i.e., X' lies to the 
right of X. We will refer to minimum- weight perfect matchings on X and X' as 
partial matchings and to the minimum-weight perfect matching on X U X' as joint 
matching. The following result is closely related to properties of "local matching 
indicators" introduced and studied in [31[S]. 

Theorem 5. Whenever an arc {xi,Xj) [respectively {x'^i,x'ji)] is hidden in the par- 
tial matching on X [respectively on X'], it belongs to the joint optimal matching 
and is hidden there too. 

Observe that exposed arcs in partial matchings are generally not preserved in 
the joint matching: they may disappear altogether or become hidden. 

Proof. By contradiction, assume that some of hidden arcs in the partial matching 
on X do not belong to the joint matching. Then there will be at least one exposed 
arc (xe, Xr) in the partial matching on X such that some points Xi with xi < Xi < Xr 
are connected in the joint matching to points outside 

Indeed, if points inside every exposed arc {xg, Xr) would be matched in the joint 
matching only among themselves, then by LemmaUtheir matching would be exactly 
the same as in the partial matching on X, and therefore all hidden arcs between Xi 
and Xr would be preserved in the joint matching. 

Suppose {xi, Xr) is the leftmost arc of the above kind. Denote all the points in 
the segment [xi,Xr] that are connected in the joint matching to points on the left 
of xi hy zi < Z2 < ■ ■ ■ < Zk', denote the opposite endpoints of the corresponding 
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arcs hy yi > y2 > ■ ■ ■ > yk, where the inequahties follow from the fact that the joint 
matching is nested. Likewise denote those points from [x£,Xr] that are connected 
in the joint matching to points on the right of Xr hy z[ > z'2 > ■ ■ ■ > z^, and their 
counterparts in the joint matching hy y[ < y2 < ■ ■ ■ < y'f., (fig- HI). 

Although k or k' may be zero, the number k + k' must be positive and even. 
Indeed, by Corollary [3] the segment [x^jX^] contains an even number of points 
and all of them must be matched in a perfect matching; removing from the joint 
matching all arcs whose ends both lie in [xe,Xr], we are left with an even number 
of points that are matched outside this segment. 

Let us now restrict our attention to the segment [xg, Xr] and consider a matching 
that consists of the following arcs: those arcs of the joint matching whose both 
ends belong to [a;£,x,-]; the arcs (zi,Z2), {z2k-i, Z2k), whercQ k = [A;/2J; the 
arcs (z2j . . . , (zg^,, Z2^t_i), where k' — [k' /2\ ; and (zfe, z^,) if both k and k' are 
odd. 

Denote by W' the weight of this matching. By Lemma IH it cannot be smaller 
than the weight Wq of the restriction of the partial matching on X to [a;^,^^]. For 
the total weight W of the joint matching on X U X' we thus have 

(4) W >W -W + Wo. 

The right-hand side of (|4]) is represented in fig. [2] (a) in the case when both k 
and k' are odd. It is a sum of positive terms corresponding to the arcs of the joint 
matching outside [xg, Xr] (not shown), the arcs of the partial matching on X inside 
[a;^,!^] (not shown, with exception of {x£,Xr) represented with a solid arc in the 
upper halfplane), the arcs of the joint matching having one end inside [xg, Xr] and 
the other end outside this segment (solid arcs in the upper halfplane), and negative 
terms that come from subtraction of W and correspond to the arcs connecting the 
z points (solid arcs in the lower halfplane). 

We now show that by a suitable sequence of uncrossings the right-hand side 
of ([4|) can be further reduced to a matching whose weight is strictly less than W. 

Step 1. Note that the arcs (zi,?/i) and {xi,Xr), shown in fig. [2] (a) with thick 
lines, are crossing. Therefore 

d{yi,zi) + d{xi,Xr) > d[yi,xe) + d{zi,Xr). 

Uncrossing these arcs gives the matching represented in fig. [2] (b) and strictly re- 
duces the right-hand side of (H]): 

W >W - W' + Wli- %i, zi) - d{xt,Xr) + d{yi,xi) + d{zi,Xr). 

Now the arcs (2/2,^2) and {zi,Xr) are crossing, so 

d{y2,Z2) + d{zi,Xr) - d{zi,Z2) > d{y2,Xr) 

and therefore 

W > W-W' + W[,-d{yi,zi)-d{y2,Z2)-d{xe,Xr) + d{yi,xe) + d{zi,Z2) + d{y2,Xr)- 
The right-hand side of this inequality is represented in fig. [2] (c). 



[^J is the largest integer n such that n < ^. 
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Vk ■■■ yi xe 



Zk ■ • ■ Zi Xy yi 



(c) 



Figure 2. Step 1 of the proof (see explanation in the text). 



Repeating this step k = [k/2\ times gives the inequahty 



W>W-W' + W^-d{xe,Xr)- YI diVi'^i) 

l<i<2K 

l<i<K l<i<K 



where in the rightmost sum yo is defined to be xe. Note that at this stage all arcs 
coming to points zi,Z2,--- from outside [x^, Xr] are eliminated from the matching, 
except possibly {yk,Zk) if k is odd. 

Step 2. It is now clear by symmetry that a similar reduction step can be 
performed on arcs going from z'i,z'2, ... to the right. Repeating this k' = [k'/2\ 
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(c) 

Figure 3. Step 3 of the proof. Note that in stage (c) the arc 
(zfcjZ^,) gives two contributions with positive and negative signs, 
which cancel out each other. 



times gives the inequaUty 

W>W -W + w;,-d{xt.xr)- d[y^,z,)- M',y[') 

l<i<2K l<i'<2K' 
+ Yj d,{z2i-l,Z2t) + Y d{y2i~i,y2i-2) 

+ Y '^(4^'>4i'-l) + Y '^iy2i'-2^y2i'-l)+diy2n,y2n') 
l<'i'<^' l<i'<K' 

where y'g = Xr- 

Step 3. If k and k' are odd, we perform two more uncrossings shown in fig. 
The final estimate for W has the form 

(5) W>W-W' + W^-d{xi,Xr)- J2 E ^^(4,^^') 

l<i<k l<i'<k' 

+ E d{z2i-i, Z2i) + E d{z2i, , z'2^,_i) + d{zk, z'^.,) ■ [k, k' are odd] 
+ Y d{y2i-i,y2i-2) + Y d.iy2t'-2,y2t'~i) + d{yk,y'k') ■ [k, k' aie even] 

l<i<.K l^i'^^' 

where notation such as [A;, k' are odd] means 1 if A;, A;' are odd and otherwise 
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The right-hand side of ([5]) contains four groups of terms: first, 
W - d{y^,z,)- d{z[,,y[,), 

l<i<k l<i'<k' 

corresponding to the joint matching without the arcs connecting points inside 
[a;£,a;r] to points outside this segment; second, 

W- Y d.{z2i-i,Z2i) - Y d{z'2,,,Z2i,_-^) ~ d{zk,zl,) ■ [k, k' a.Te odd], 

l<z<At l<i' <k' 

which comes with a negative sign and corresponds to the arcs of the joint matching 
with both ends inside [xe, Xr], and cancels them from the total; third, 

Wq - d{xi,Xr), 

with positive sign, which corresponds to the hidden arcs of the partial match- 
ing on X inside the exposed arc not including the latter; and finally the 
terms in the last line of ([5]), corresponding to the arcs matching xe, a;,., and points 
Ui, . . . ,yk,y'i, ■ ■ ■ ,y'i^, , i.e., those points outside [xi,Xr] that were connected in the 
joint matching to points inside this segment. 

Gathering together contributions of these four groups of terms, we observe that 
all negative terms cancel out and what is left corresponds to a perfect matching 
with a weight strictly smaller than W, in which all arcs hidden by {xe,Xr) in the 
partial matching on X are restored. There may still be some crossings caused by 
terms of the fourth group and not involving the hidden arcs in [a;^,^:^]; uncrossing 
them if necessary gives a nested perfect matching whose weight is strictly less than 
that of the joint matching. This contradicts the assumption that the latter is the 
minimum- weight matching on X U X' . Therefore all hidden arcs in the partial 
matching on X (and, by symmetry, those in the partial matching on X') belong to 
the joint matching. □ 

5. Recursion for minimum weights 

In this section we show how to apply Theorem [5] to compute the minimum- 
weight perfect matching algorithmically. Let xi < X2 < ■ ■ ■ < X2n be a set of points 
on the real line R equipped with a homogeneous distance d of the form ([T]). For 
indices i, j of opposite parity and such that i < j, let Wij be the weight of the 
minimum- weight perfect matching on the j — i + I points Xi < x^+i < ■ ■ ■ < Xj. It 
is convenient to organize weights Wij with i < j into a pyramidal table: 

Wi,2n 



Wi^e W2J VK3,8 W4,9 
WKl,4 W2,5 VK3,6 VK4,7 W^5,8 
Wi^2 W2,3 W3,4 VK4,5 Ws^e VKfiJ 

Theorem 6. For all indices i, j of opposite parity with 1 < i < j < 2n, weights Wij 
satisfy the following second-order recursion 

(6) Wij = min [d{xi,Xj) + Wi+ij-i, Wij-2 + Wi+2,j - Wi+2,j-2] 

with "initial conditions" 



W2n-3,2n 



(7) 



Wi,i_i = 0, Wi+2,j-i = ~d{xi,Xi+i). 
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Proof. By an abuse of notation, we will refer to the minimum- weight perfect match- 
ing on points Xr < Xr+i < ■ ■ • < Xs slS the "matching Wr,s-" 

Consider first the matching that consists of the arc {xi , Xj ) and all arcs of the 
matching M^i+i.j-i, and observe that by Lemma|4]its weight d{xi,Xj) + Wi+i_j-i 
is minimal among all matchings that contain (xi ,Xj). 

We now examine the meaning of the expression + Wi+2,j — W^i+2.j-2- 

Denote the point connected in the matching Wi,j-2 to Xi by Xk and the point 
connected to x^+i by xg. By Corollary [31 the pairs of indices z, k and i l^H both 
have opposite parity. Assume first that 

(8) Xi+i < Xi < Xk < Xj-2- 

Applying Theorem [5] to the sets X = {xi,Xi+i} and X' = {2:^+2, and 
taking into account parity of k and we see that Xk and xe (as well as their 
neighbors Xk+i and xe-i if they are contained in [xi+2, Xj-2]) belong to exposed arcs 
of the matching Wi+2.j-2- Thus the matching Wi,j-2 has the following structure: 




where dashed (resp., dotted) arcs correspond to those exposed arcs of the match- 
ing Wi+2,j-2 that belong (resp., do not belong) to Wij-2- 

Since points x^-i and Xk+i belong to exposed arcs in the matching Wi+2.j-2, by 
Lemma[4]we see that the (possibly empty) parts of this matching that correspond to 
points Xi^2 < ■ ■ ■ < xi-i and Xk+i < ■ ■ ■ < Xj^2 coincide with the (possibly empty) 
matchings Wi+2,1-1 and Wk+i,j-2- For the same reason the (possibly empty) part 
of the matching Wij-2 supported on Xf+i < • • ■ < Xk-i coincides with We+i^k-i- 
Therefore 

(9) Wij-2 = d{xi,Xk) + d{xi+i,x() + Wi+2,e-i + We+i^k-i + Wfe+i,j-2- 

Taking into account ([7]), observe that in the case k = i + \ and £ — i, which was 
left out in this expression still gives the correct formula Wij-2 = d{xi,Xi+i) + 

Wi+2j-2- 

Now assume that in the matching Wi^ij the point Xj is connected to Xf and 
the point Xj-i to Xk'- A similar argument gives 

(10) W.i+2,j = Wi+2,<;'-i + Wt+i.k'-i + Wk'+ij-2 + d{xi,,Xj) + d{xk' , xj^i); 

in particular, \i = j — \ and k' = j, then Wi+2,j = W^i+2.j-2 + d(a;j-i, x^). 

Suppose that Xk < Xi>. Using again Lemma [4] and taking into account that Xk, 
Xk+i, xt~i, and xpj all belong to exposed arcs in M^i+2.j-2, we can write 

(11) Wk+l,j-2 Wk+l,l'-l + VF"rj-2, W^+2.t'-l Wi+2,k + Wk+ll'-l 

and 

(12) W^+2,0-2 = Wi+2,k + Wk+lJJ-l + VK£'j_2. 
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Substituting ([TT|) into (|9]) and (fTO|) and taking into account (fT2|) . wc obtain 

+ Wk+i,i'-i + d{xi>,Xj) + Wt'+i.k'-i + d{xk',x.j-i) + Wk'+i.j-2- 

The right-hand side of this expression corresponds to a matching that coincides with 
Wij_2 on [xi,Xk], with Wi+2j-2 on [xk+i,xe'~i], and with W^+ij on [xi',Xj]. By 
Lemma |4] this matching cannot be improved on any of these three segments and 
is therefore optimal among ah matchings in which Xi and Xj belong to different 
exposed arcs. 

It follows that under the assumption that Xk < xt the expression in the right- 
hand side of Q gives the minimum weight of all matchings on Xi < Xi^^i < ■ • • < 
Xj. Moreover, the only possible candidates for the optimal matching are those 
constructed above: one that corresponds to d{xi, Xj) + Wi+ij-i and one given by 
the right-hand side of the latter formula. 

It remains to consider the case Xk > xi' . Since Xk 7^ xii for parity reasons, it fol- 
lows that Xk > Xi' ; now a construction similar to the above yields a matching which 
corresponds to Wij-2+Wi+2,j^Wi+2.j-2 and in which the arcs (xi, Xk) and (xi'^Xj) 
are crossed. Uncrossing them leads to a matching with strictly smaller weight, which 
contains the arc {xi,Xj) and therefore cannot be better than d{xi,Xj) + M^i+i,j-i. 
This means that ([6]) holds in this case too with Wi^j = d{xi,Xj) + Wi+ij-i. □ 

Obviously, recursion ([6]) can be solved for all 1 < z < j < 2n in 0{n^) operations, 
resulting in computation of weights Wij of all partial optimal matchings. This 
process is carried out in a "bottom to top" fashion: in the pyramid, weights Wi,j 
with smaller values oi j — i are computed first. 

To determine the optimal matching on all the points xi,X2t ■ ■ , X2n, one should 
keep track of those pairs for which minimum in ^ is attained at the first 

alternative. Indeed, if the minimum is never attained at the first alternative, then it 
is easy to see that the optimal perfect matching is (xi,X2), (xs, 0:4), . . . , (a;2n-i, X2n)- 
Suppose now that the first alternative provides minimum for for some Wi^j^ . Then 
according to Theorem[5]one can retain the matching for the points Xi„^i, . . . , Xj^^i 
that has been computed by this moment and consider a new, smaller minimum- 
weight perfect matching problem on the points xi , X2 , . . . , , , , . . . , X2n ■ 
More precisely, it suffices to remove from the pyramidal table quantities Wij with 
i,j satisfying at least one of the conditions ig < i < jo or io < j < jo, replace Wi^j^ 
with d{xig, Xjg), stack the cells with either i or j outside {io,jo) into a smaller 
pyramidal table, and continue solving the recursion. 

This reduction step is illustrated in fig. HI The elements to be removed from the 
table are shown in light gray in the top pane. Assuming that the rows are scanned 
left to right, at the moment when it is found that the element Wi^^j^^ = W^^ (in 
brackets) involves the first alternative in (|6]), the elements shown in parentheses 
have not yet been computed. Those of them that are to be kept in the table (the 
"black" ones) are then stuck with the already computed elements to form a smaller 
pyramid shown in the bottom, and the recursion resumes with the element 
replaced with d{x3,xe). 

At the reduction step illustrated in fig. |4l the matching is updated with one 
arc {x4, xs), which is guaranteed to belong to the optimal matching by Theorem [5l 
Generally, every time a reduction step is performed on an element Wi^jg and results 
in removal of the points Xk^ < xu^ < ■ ■ ■ < a;fe2„ that have been "covered" with 



12 



J. DELON, J. SALOMON, AND A. SOBOLEVSKI 



(VKi.io) 

(W^i,8) (T^2,9) (Wa.io) 

(VFl,6) (W^2,7) (W^3,8) iVn.gj (Ws.io) 

[^-3,6] (VK5.8) (VFe.g) (W^7,io) 

W^l,2 M^2,3 W^3,4 ^4,5 W5,6 We,r W7,8 W8,9 Wg.io 

i 

(Wi.s) {W2.9) {W3.10) 

(W^i.e) {W2,7) {W3,s) (W^e.g) (l^rao) 

W^i,2 W^2,3 rf(a:3,a;6) 14^6,7 Wjs Ws,9 Wg,io 

Figure 4. Reduction of the pyramidal table (iq = 3, jo = 6, 71 = 5). 

the arc (io,jo), the matching is updated with arcs {xk-i^,Xk2), ■ ■ ■ A^k2m-n^k2m)- 
When the process reaches the top element of the pyramid, a similar matching of 
consecutive pairs is performed on the remaining elements. 

A formal description of this procedure is given in in somewhat different 

terms (emphasizing "local matching indicators" instead of the quantities Wij). 
Since this algorithm requires computing all the weights Wij in the worst case, 
it takes at most 0{n^) operations for the matching on 2n points. 

It is interesting to compare performance of our algorithm to the celebrated Ed- 
monds algorithm for minimum- weight perfect matching in a general graph with N 
vertices and M edges [BJ[7|. Note that we effectively consider a complete graph 
on = 2ri vertices, i.e., M = n{2n — 1). In this case recent improvements of 
the Edmonds algorithm run in 0{n^) time (see, e.g., the survey [3]) whereas our 
specialized algorithm is at most quadratic in n. 
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